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THREE-DIMENSIONAL  SIMULATION  OF  STRESS 
AROUND  INCLUSIONS  IN 
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{Receired  15  biovemher  1989.  and  in  hnal  form  5  June  1992) 


A  (initc  element  numerical  model  is  presented  which  simulates  vibrations  in  and  around 
inclusions  in  vi.scoelastic  tluids  excited  by  a  sinusoidal  pressure  wave  or  pulse  generated  by 
noisy  machinery.  Stress  transler  processes  are  investigated  and  regions  ol  highest  energy 
dissipation  identified.  The  model  simulates  the  time  evolution  of  stress  generation  in  and 
around  inclusions  in  the  viscoelastic  plastic  fluid  in  three  space  dimensions  Fluid 
composites  recommended  (or  their  power  dissipative  characteristics  are  ranked  according 
to  their  ability  to  dissipate  acoustic  power  in  a  simple  composite  structure. 


1  INTRODUCTION 

Many  layers  of  material,  which  cover  engineering  structures  such  as  the  human  body  of 
ancchoic  and  noise  and  vibration  decoupling  coatings,  are  inhomogeneous  in  structure. 
The  inhomogeneitics  may  be  hair  roots  in  skin,  or  air  cavities  in  anechoic  layers,  which 
can  maximize  dissipation  of  vibrational  or  acoustic  energy  in  the  layer.  The  precise  manner 
tn  which  the  inhomogeneities,  particularly  the  inserts  of  air  often  found  in  anechoic 
materials,  scatter  vibration  and  acoustic  energy  is  not  well  understood.  Most  theoretical 
approaches  w  hich  have  addressed  these  problems  to  date  have  used  fairly  simple  analytical 
approaches  to  describe  the  mechanisms  thought  to  be  involved  [1).  The  present  study  aims 
to  complement  the  analytical  work  by  developing  a  numerical  scheme  which  may  provide 
a  technique  for  elucidating  the  areas  of  uncertainity  in  the  analytical  studies  of  the  acoustic 
response  of  materials  with  insets  and  also  media  excited  by  an  insert,  such  as  skin  pricked 
by  a  hair. 

A  primitive  finite  diflcrence  model  is  developed  which  is  able  to  simulate  the  vibration 
and  acoustic  response  of  a  fluid  to  an  incident  pressure  wave.  The  fluid  may  have  elastic, 
viscous  and  plastic  properties  and  contain  inclusions  which  may  be  air  or  fluid. 

The  fluid  is  modelled  as  consisting  of  a  large  number  of  cubical  elements  which  may  be 
distorted  by  the  progression  of  an  acoustic  or  vibration  excitation.  The  model  is  able  to 
simulate  a  limited  degree  of  macroscopic  distortion  and  able  to  trace  the  passage  of  a 
pressure  pulse  or  sinusoidal  wave  through  a  multi-layer,  multi-cavity  or  multi-inclusion 
fluid,  such  as  coatings  used  to  quieten  metal  structures  or  an  elastomeric  mount  designed 
to  attenuate  the  transmission  of  sound  and  vibrations  over  a  wide  frequency  range, 
including  frequencies  at  which  sections  of  the  mount  resonate  [2].  The  results  of 
representative  calculations  show  the  response  of  cubical  cavities  and  inserts  embedded  in 
viscoelastic  or  plastic  media  to  a  sinusoidal  excitation.  The  acoustic  response  of  cylindrical 
cavities  in  an  elastic  fluid  is  also  shown.  The  model  has  also  been  used  to  simulate  the 
response  of  surrounds  of  double  cone,  diamond  and  ellipsoidal  cavities. 
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Figure  I .  Section  of  rectangular  column  cavity  modelled  and  typical  surround  of  cavities.  The  unit  cell  is  shown 
by 


As  many  fluids  used  to  quieten  metal  structures  are  backed  by  steel  plates,  the  model 
considers  the  clastic  behaviour  of  a  steel  backing  plate  in  calculating  the  acoustic  or 
vibrational  response  properties  of  the  fluid  systems  described.  The  model  permits 
consideration  of  water  in  front  of  the  fluid  and  air  or  water  behind  the  steel  backing  plate. 

Typical  cross-sections  of  the  fluid  coatings  which  have  been  modelled  are  shown  in 
Figure  1  and  types  of  inserts  and  their  pattern  in  Figure  2.  The  model  is  general,  however, 
and  is  able  to  represent  a  wide  variety  of  acoustical  and  vibrational  responses  of  materials 
which  may  be  important  for  anechoic  coatings,  anti-vibration  and  acoustic  isolation 
mounts  for  machinery  and  the  prickle  response  of  the  human  skin. 


2.  MODEL 


2.1.  PROPAGATION  MODEL 

The  propagation  of  a  pressure  or  stress  wave  through  a  compressible  medium,  such  as 
the  fluids  considered,  is  described  by  the  laws  of  conservation  of  mass  and  momentum  and 
the  equation  of  state  which  relates  the  pressure  or  stress  in  the  medium  to  the  strain  and 
its  material  properties  [3]. 

The  displacement  of  fluid  particles  by  an  acoustic  wave  is  small  relative  to  the 
characteristic  dimensions  of  a  typical  acoustic  structure  such  as  an  anechoic  coating  on 
a  steel  plate.  The  behaviour  of  such  relatively  weakly  perturbed  systems  can  be  solved  by 
numerical  methods  by  using  the  Lagrangian  technique  [4].  This  technique  is  particularly 
suited  to  calculations  in  which  discontinuities  in  physical  properties,  such  as  density  and 
velocity,  occur. 

2.2.  CONSERVATION  OF  MASS 

The  law  of  conservation  of  mass  for  a  compressible  medium  is  usually  expressed  in  an 
Eulerian  framework  as  follows:  the  time  rate  change  of  mass  density  at  any  point  is  equal 
to  the  negative  divergence  of  the  momentum  at  that  point.  This  law  can  be  expressed  in 
mathematical  form  in  a  frame  of  reference  which  moves  with  the  fluid,  the  Lagrangian 
frame,  and  is 


6p  dt  =  —  pV  •  V. 


(1) 
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Figure  2  Arrays  of  typical  inserts  or  cavities.  Various  unit  cells  are  shown  by - and  labelled  a.  b  and  c. 

respectively. 


where  p  is  the  mass  density  and  v  is  the  local  velocity  of  the  fluid.  The  discrete  numeri¬ 
cal  representation  of  the  Lagrangian  approach  defines  cells  of  material  corners,  and  hence 
boundaries  which  move  with  the  local  fluid  velocity.  In  the  present  study  three- 
dimensional  (.V.  y. : )  space  is  considered  to  examine  details  of  shear  and  normal  stresses 
in  the  three  orthogonal  directions.  In  this  system  the  indices  i,  /  and  k  indicate  cell  position 
counters  in  the  .v,  »■  and  r  directions,  respectively.  To  consider  firstly  motion  in  the  .v 
direction,  the  position  of  a  cell  corner  at  time  i  +  di.  -i-  is  related  to  its  previous 
position  .v„*(0  at  a  time  t  by 


x,^{i  +  St  )  =  x,^(t)  + 


t*  fit 


(2) 


If  the  velocity  variation  over  the  interval  St  is  represented  by  an  average,  the  above 
expression  can  be  written  as 


I  ^  J.  ,.n  +  I 


•v;,:  =  x'u  +  i-z 


(3) 


where  the  integer  superscript  n  indicates  a  particular  instant  of  time  via  the  relationship 
t  =  n  St. 

Corresponding  expressions  can  be  derived  for  the  y  and  z  co-ordinates  of  the  cell  corner. 
In  the  time  interval  —  I  to  n.  the  density  variation  of  the  cell  designated  ijk  which  is 
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Figure  Vectors  of  dement  of  volume 


bounded  by  the  corner  indices  i/k.  i  -f-  I/A.  i  +  1/  +  U',  i  +  I./  +  lA  +  I.  //  +  lA.  //  +  lA  4-  I, 
ijk  +  1,1  +  1;A  +  1,  is  upon  ustnc  a  vector  cross-product  notation,  given  by 

/'  =  I  ^  b)|  +  |f-(d  X  e)l!"  ' 

^  ■  !lc  ■  (a  X  b)l  +  If  •  (d  X  e)!!" 

where  a,  b.  c.  d,  c  and  f  are  the  vectors  of  the  side  of  the  cell  shown  in  Figure  3.  and  where, 
for  example,  a  is  the  vector  v.  i,  2|„,  and  .v,  r,  r|,, 


2.3.  CONSERVATION  (»l  MOMI  NTI  M 

The  law  of  conservation  of  momentum  is  commonly  expressed  in  the  moving, 
Lagrangian  frame  as 

/)  dv  d/  =  —  Vp.  (5) 

In  a  solid  material  the  stress,  rather  than  pressure,  causes  the  cell  boundaries  to  move  and 
rotate  (5).  and  the  conservation  of  momentum  is  expressed  in  terms  of  the  stress  a,j  and 
the  velocity  components  along  the  three  axes  of  Cartesian  co-ordinates  in  tensor  notation 
[6]  as 

P  d/  ,  d/ =cff„/c.v,.  (6) 


By  extending  the  one-dimensional  situation  defined  by  equation  (2)  to  the  three-dimen¬ 
sional  case,  the  velocities  of  the  cell  corners  are  obtained  from  equation  (6).  with  the 
discretization  of  equation  (3).  These  arc.  for  i\ 
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and  for  i . , , 


J,<<*  I  :  _  -  I 

-/4  •i/i 
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III  pL. 


'  >  a'l,  —  a., 

♦••Hi  —4*1*  -  1 

/  m  i  -  I  m  •  /  -  I  ^  tmk  ^  imk  -  I 


+  i 


k 


I 


*  l/f  ^  h  ~  lo 


+  I 


(9) 


Here  )  is  the  position  vector  of  the  centre  of  the  cell  designated  ijk  at  time 

increment  /i,  the  components  of  which  are  defined  by 


/-flr-lil-fl 

i  I  ^"ik~iY.  1 1 I  i 

I m  I  m  ^  f  n  »  k  Irntm^fomk  I m  I  m  m  f  o  m  k 

This  system  of  equations  is  complemented  by  the  equation  of  state  of  the  material,  which 
links  the  strain  experienced  by  the  fluid  of  the  stresses  imposed  on  it,  to  provide  a  complete 
description  of  the  response  of  a  material  to  an  acoustical  or  vibrational  excitation. 


2.4.  STRAIN  DISPLACEMENT  RELATIONS 

The  strains  are  determined  from  the  displacement  of  a  particle  in  a  solid  or  fluid  [6]. 
It  is  convenient  to  consider  the  particle  at  position  .y"*  at  time  t  of  index  n  to  have  already 
received  a  small  displacement  u"  ^  in  the  v  direction  so  that  its  undisturbed  position  was 
at  Then,  in  the  present  finite  primitive  element  model  the  displacements  along 

the  three  axes  of  the  Cartesian  co-ordinate  system  are  defined  by 


Sm  ''//A  Ifk  ^  }$/k*  •ttk' 


(ID 


The  incremental  change  of  these  displacements  leads  to  the  deformation  of  an  element  of 
volume  and  can  be  expressed  mathematically  in  terms  of  the  divergence  and  curl  [10]  of 
the  displacement.  For  instance,  the  deformation  of  the  element  of  volume  in  the  .y  direction 
can  be  expressed,  by  the  chain  rule  [6,  7],  as 


(12) 


This  may  be  expressed  in  mathematical  shorthand  by 


(13) 


where  higher  order  terms  in  dx"^,  and  have  been  neglected.  The  and 
represent  rotation.  The  and  represent  the  direct  and  symmetric  strains 

needed  to  compute  the  stresses  induced  in  the  material.  The  strains  at  the  cell  centres  of 
the  number  scheme  used  here  are  related  to  the  displacements  at  the  cell  corners  by 


-u’i 

,  _  I  ^  y  *■  ♦ 

‘  “  4  iL  Z,  n  _ 

m  •  f  o  ^  k  I  It 


1  i  +  I  y"  —  W; 


c"  =1  y  y  -III 


O+Io  -'■//<> 


i  -I  M?  —u" 

("  =  i  y  y  •*"*  * ' 


(14-16) 
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'  - 1  ‘  '  m"  -  w  I  - 1  *  - 1  i/"  -  w:  \ 

k  •'  1  \nu>  i/fut  /»//**  i  -r  |««  .*  //.•  / 

(17) 

-W?  u"  -u"  \ 

I  I  -7^  — +  I  I 

#  »i  •  4'  1  ♦  lnt4>  *'  mill  1  =  i  m  -  >  •  tmk  ♦  I  Imk  / 

(1«) 

'•  1  ‘  '  1/?  -  -  -  1  1  u"  -  u"  \ 

/  A  1  n  »  A  J  //  f  hi  /ft!  imi  -•  1  •  /w4  / 

(19) 

The  displacemenis  at  cell  corners  are  calculated  by  integrating  the  velocities  at  the  corners 
via  equation  (3). 

2.5.  STRESS- STRAIN  RELATIONS  AS  EQUATIONS  OF  STATE 

The  simple  theory  of  elasticity  assumes  that  the  induced  stresses  and  strains  are  linearly 
related  to  each  other  and  that  the  stress  and  strain  tensors  always  have  the  same  axes.  If 
the  stresses  and  strains  along  three  orthogonal  axes  are  considered,  then  the  general 
stress-  strain  relation  can  be  written  as  [6], 

ff„  =  (20) 

The  parameter  represents  the  sum  of  the  orthogonal  strains:  /  and  are  known  as 
Lame's  constants  and  //  is  also  called  the  modulus  of  rigidity  and  measures  the  resistance 
of  the  substance  to  distortions.  These  constants  are  related  to  Young's  modulus  E  and  the 
Poisson  ratio  v  by 


;.  =vJE' (!  +  »•)( I -2r).  /i=£;2(l+v).  (21,22) 

The  stress  situation  is  described  as  hydrostatic  when  the  direct  stresses  in  three  orthogonal 
directions  are  equal  and  the  shear  stresses  are  zero.  The  direct  strains  will  then  be  equal. 
Under  these  conditions,  in  general 

a,,  =  o..  =  =  ^ki..  (23) 

where  k  =/  +2/l3.  The  parameter  k  is  called  the  bulk  modulus  because  the  relative 
change  of  volume  if  3< ,,  to  first  order,  so  that  k  is  the  ratio  of  the  symmetrical  stress  to 
the  change  in  volume.  For  an  ideal  fluid,  /i  =  0  and  k  = 

In  materials  which  behave  in  a  hystereiic,  lossy  fashion  a  phase  delay  is  induced  in 
acoustic  signals  propagating  through  them.  Accordingly  such  materials  are  characterized 
by  Lame  constants  involving  differential  operators  with  respect  to  time  in  the  time  domain 
and  complex  constants  in  the  frequency  domain  [I]. 

Equation  (20)  decomposes  in  the  orthogonal  co-ordinate  (x,y.z)  system  used  in  the 
present  simulation  model  at  any  instant  of  time  defined  by  the  time  index  n  to 

ff:,  =  (/  -(-2(i)(':,-i-/.(r:,-i-f?,)  (24) 

and 

a",,  =  u  +  2g)(';.  -f  -I-  (^,),  (T';,  =  (;.  +  2n)c".,  -t-  /M:.  +  r :,  )■  (25.  26) 

Also,  a",  a",  =  2n(”, .  a",  =  <t”.  =  2/i(",  and  it",  =  cr".  =  2/i(“.,  where  /  and  (x  involve 

differential  operators  with  respect  to  time  for  lossy  materials. 
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2.(1.  STRESS  STRAIN  RELATIONS  i  OR  LINEAR  PLASTIC  -  VISCOELASTIC  MATERIALS 

The  vibrational  behaviour  ol  many  fluids,  or  materials  which  behave  like  fluids  such  as 
elastomers,  varies  with  temperature,  pressure  and  rate  of  strain.  The  elastic  behaviour  may 
vary  when  stresses  are  repeatedly  reversed  and  the  material  may  exhibit  hysteresis.  Once 
the  yield  point  of  a  substance  exhibiting  plastic  behaviour  is  passed,  increasing  stresses  give 
steadily  increasing  strains.  Hence  any  single  mathematical  model  can  be  expected  only  to 
approximate  the  elastic  behaviour  of  actual  substances  under  limited  conditions  [5].  The 
behaviour  of  all  the  fluids  considered  in  this  study  is.  to  a  good  approximation,  linear  for 
the  small  pressure  fluctuations  considered  to  be  typical  of  acoustic  and  vibration  generated 
by  machinery. 

The  simplest  general  linear  relationship  between  a  change  in  stress  and  associated  strain 
can  be  written  in  tensor  notation  as 


fT,  +  /„  rn,i  ri  =/.(!+  t^  r  rt  +  2/i(  1  +  /;.  r  i  t )(.  .  (27) 

where  /,  and  t.  are  characteristic  relaxation  times  of  the  material. 

Plastic  behaviour  is  described  by  the  condition  that,  for  a  constant  stress  a,,  applied  at 
t  =  0.  the  material  response  can  be  characterized  as  follows  [8]; 

fT,,  +  ((T,  -  a„)t  tx- /.{\  -E cr  )<„,„,<>„+ 2//(  I +/:  r/c'/ )«,,  if(j,,^or„.  (28) 

The  /:  and  /,  are  characteristic  relaxation  times  of  the  material,  and  describe  the  complex 
nature  of  the  clastic  moduli  in  the  time  domain  used  commonly  in  steady  state  impedance 
fomulations. 

The  linear  model  of  equations  (27)  and  (28)  gives  a  representation  of  damping  of 
vibrations  by  internal  friction.  The  internal  friction  causes  a  phase  delay  in  the  trans¬ 
mission  of  steady  sinusoidal  signals  through  the  material  which  can  be  expressed  as  a 
loss-tangent,  tan  (').  which  is  related  to  the  relaxation  times  and  the  frequency  lo  of  the 
signal  by 


tan  <),=  w(/| —  /„),'(  1  4- oj'/,  r,,).  tan  =  ty((;  — /o)/(  1  +  wV^o)-  (29.30) 

The  angle  6  measures  the  lag  of  strain  behind  stress,  and  is  known  as  the  loss  angle  of 
the  material  and  provides  a  measure  of  the  internal  damping  of  stress  waves  [5]. 

The  loss  tangents  of  materials  commonly  used  in  acoustic  structures  cannot  always  be 
related  to  the  excitation  frequency  as  described  by  equations  (29)  and  (.^0).  The  magnitudes 
of  /,  and  N  vary  with  frequency  such  that  a  simpler  Kelvin  Voigt  model  [5).  where  /„  =  0. 
might  be  adequate  to  describe  the  ancchoic  coating  under  consideration.  The  more  general 
formalism  of  equations  (27)  and  (28)  is  kept  to  widen  the  range  of  acoustic  structures  which 
can  be  studied  by  the  present  model.  In  particular,  many  elastomers  are  characterized  by 
different  values  for  the  loss  tangent  associated  with  the  first  and  second  Lame  constants 
which  is  reflected  in  the  time  differential  operator  coefficients  of  these  constants  in 
equations  (27)  and  (28). 

The  viscoelastic  stress-strain  equation  can  be  expressed  in  a  finite  cell  formulation  which 
relates  the  stress  tensor  a,,  at  time  index  n  and  cell  centre  {i  jk )  to  the  corresponding  strain 
tensor  arising  from  the  movement  of  the  adjoining  cell  corners  (see  Figure  4)  by  using 
backward  differences  for  the  time  step,  where  the  <  s  are  calculated  from  the  movement 
at  cell  corners.  The  numerical  procedure  for  evaluating  velocity  and  displacement  at  the 
cell  corners  and  density  and  stresses  and  strain  at  cell  centres  enhances  the  computational 
stability  of  the  Lagrangian  technique  used. 
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2.7,  I  OSS  I  rN(  iioN 

Power  dissipation  considered  is  caused  by  loss  mechanisms  which  become  active  during 
the  delormation  ol  ihe  material.  The  power  dissipated  b\  these  mechanisms  is  given  b> 
dH  d/,  where  H  is  the  work  done  by  the  stresses  across  the  surface  of  an  element  of 
volume  ()T  =  <)v<>r<)c.  and  the  forces  /»rr,  rt.  fu'i,  ri  and  pfr,  ci. 

For  an  infinitesimal  increment  of  displacement  the  work  increment  per  unit  volume  is 
given,  in  tensor  notation  for  simplicity  [3].  by 

=  .  (31) 

where  the  summation  signs  may  be  e>mtltcd  under  the  summation  convention  for  tensors 
[9].  The  work  increment  hM  and  the  time  increment  <)i  can  be  related  by 

()(„  <>/.  (32) 

where  the  instantaneous  power  absorbed  is  given  by  dP'/dr  as  Pi  approaches  zero. 

U  the  incident  wave  is  steady  state  cyclic  with  angular  frequency  tu.  the  average  power 
dissipated  per  unit  volume  is  of  interest.  This  can  be  written  as 


(.1  r'--’-dM- 

—  --dr 
2ndr  J,  di 


Taking  limits  in  equations  (32)  and  substituting  m  equation  (33)  gives 

L  =  f'  "  a,  dr  (34) 

2?:  !,  dr 

*  ' 

where  /  is  the  time  at  the  beginning  of  the  cycle,  t  +  2n  v)  is  the  time  at  the  end  of  the 
cycle  and  2n  id  is  the  period  of  the  cycle.  In  the  numerical  scheme  used  in  the  present  work, 
this  loss  function  L  at  all  cell  locations  ijk  is  calculated  from 

L  V  a".  (("  -i"  '  ). 


The  units  of  L  is  watt  per  cubic  metre  (W  m'i. 


(35) 
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The  model  was  programmed  to  simulate  the  acoustic  and  vibrational  response  of  fluid 
layers  containing  several  types  of  cavities  or  inserts  to  investigate  the  effect  the  cavity  or 
insert  shape  has  on  the  global  acoustic  properties  of  such  layers.  The  mam  shapes  of 
cavities  or  inserts  investigated  were  columns,  cubes  and  short  cylinders. 

The  cavities  or  inserts  were  considered  to  be  spaced  at  regular  intervals  which  permits 
the  consideration  of  a  unit  cell  only.  Typical  unit  cells  containing  a  portion  or  all  of  a  cavity 
or  insert  are  shown  in  Figures  1  and  2  by  dashes.  The  cavities  or  inserts  were  5  mm  in 
diameter  or  square  and  spaced  5  mm  apart.  The  acoustic  structure  wa.^  assumed  to  be  made 
up  of  a  layer  of  such  unit  cells  extending  to  infinity  in  the  r  and  r  direction.  Each  unit 
cell  w  as  assumed  to  behave  in  an  identical  manner  w  ith  continuity  of  excitation  across  unit 
cell  interfaces. 

The  calculations  were  carried  out  in  the  order  shown  in  the  flowchart  of  Figure  5.  A 
sinusoidal  excitation  in  the  v  component  of  the  \elocity  vector  was  started  at  plane  i  =  I 
at  time  /  =  0  or  time  increment  n  -  0.  The  displacements  of  the  cell  corners  were  calculated 
from  equation  (.^)  using  equation's  (7)  (9)  A  cross-section  of  the  primitive  finite  clement 
grid  used  in  the  present  model  is  show  n  in  Figure  4.  The  volumes  of  the  cells  and  local 
density  of  material  was  updated.  Finally,  the  stresses  were  updated  to  complete  one  cycle 


I 

All  quuniKie.s  | 
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Over  complete  i 

field  at  particular  1 
tunc  cunsiaiii 
specified 


L 


Kvaluaied  over 
one  cycle  when 
steady  state 
condition  reached 


Tune  iteration 
continued  until 
steady  state 
condition 
achieved  (eight 
cycles) 
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of  iicraiion.  The  acoustic  structures  considered  required  about  five  cycles  of  excitation 
before  steady  state  vibrational  responses  were  observed.  Computer  memory  and  comput¬ 
ing  time  considerations  limited  the  number  of  cells  in  the  calculations  to  60  along  each 
axis  of  the  Cartesian  co-ordinate  system. 

Simulations  were  first  carried  out  for  both  loss-less  and  lossy  elastomer  half-spaces 
containing  no  cavities  or  inserts  to  check  the  accuracy  of  the  computations  against 
analytical  results.  The  tests  indicated  that  the  amplitudes  of  displacements  computed  were 
accurate  to  within  half  of  one  percent.  Calculations  were  then  carried  out  to  check  the 
extent  to  which  the  compulations  conserved  energy  in  the  system  modelled.  The  calcu¬ 
lations  were  estimated  to  conserve  energy  to  about  95  percent.  This  figure  was  considered 
good  for  the  spatial  grid  used  and  the  complexity  of  the  calculations  involved  [10].  It  was 
estimated  that  the  numerical  calculations  conserved  energy  or  power  to  about  90  percent 
for  the  system  modelled. 


4.  RESULTS 

The  acoustic  behaviour  of  sound-absorbing  or  dissipative  structures  can  conveniently 
be  described  in  terms  of  the  behaviour  of  regions  where  the  major  responses  to  the  acoustic 
excitation  occur.  These  regions  can  be  considered  to  be  elements  of  the  overall  acoustic 
structure  [Ij. 

The  details  of  the  vibrational  and  acoustic  response  of  material  around  and  in  inclusions 
is  most  readily  presented  b>  considering  first  the  two  space  dimensional  response  of 
columnar  inserts  in  a  medium,  such  as  the  rectangular  columns  shown  in  Figure  I  and 
similar  cylindrical  columns  where  the  vibrational  excitation  is  perpendicular  or  normal  to 
the  long  axis  of  the  column. 

The  acoustic  response  of  the  background  medium  is  most  clearly  presented  by 
considering  first  the  case  in  w  hich  the  speed  of  sound  in  the  insert  is  very  small  compared 
with  the  speed  of  sound  in  the  background  medium.  Such  a  composite  scatters  the  sound 
at  the  acoustically  soft  surface  of  the  insert  and  facilitates  its  absorption  though  mode 
conversion,  where  the  shear  modulus  of  the  medium  enhances  the  absorption  of  the 
acoustic  excitation.  The  insert  material  with  the  lowest  practical  speed  of  sound  is  air.  and 
hence  the  stresses  generated  around  air  inserts  or  cavities  embedded  in  a  viscoelastic 
medium  will  be  presented  first.  When  a  dilatational  wave  encounters  such  cavities,  the  wave 
is  reflected  by  the  acoustically  .soft  boundary  at  the  cavity  surface  [II]  and  some 
transformation  of  normal  stress  into  transverse  and  shear  stress  occurs. 

4.1.  inti:r-ins!:rt  rigkin  ami  insi  rt  i.ik;»;s 

4.1.1.  hilcr -inxcri  column 

The  generation  of  shear  stress  is  widely  considered  to  be  an  indication  of  acoustic 
absorption  efficiency  [12]  and  hence  computational  simulations  of  shear  around  inserts  in 
vi.scoelastic  layers  arc  shown  in  Figures  6  12.  The  acoustic  response  of  a  column  of 
rectangular  cross-section  is  illustrated  by  patterns  of  shear  stresses  shown  in  Figure  6.  The 
coaling  material  is  forced  into  each  cavity  by  the  transformation  of  normal  into  shear  stress 
through  Poisson's  effect.  This  causes  the  pattern  of  shear  stress  in  the  inter-cavity  region 
to  have  a  node  along  the  centre  of  the  column.  Symmetry  and  conservation  of  angular 
momentum  require  that  edge>  on  opposite  sides  of  the  material  column  between  the 
cavities  rotate  in  opposite  directions.  The  stresses  are  concentrated  at  the  corners  of  the 
cavities  facing  the  acoustic  excitations  in  the  examples  shown  because  mode  relaxation 
removes  energy  from  the  mcideni  excitation  before  it  reaches  the  corners  further  away;  and 
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I'lpurc  (i  Shear  stress  around  eavio  edges  of  columnar  eaviis  showing  phasing  rclalionship.  Excilalion 
ineidenl  from  lop.  waler-side  of  ligure  ..  ong  v-axis. 


viscoelasiic  losses  in  the  maicrial  in  these  examples  (see  "T^^ble  I)  have  also  caused 
attenuation  of  the  acoustic  wase  during  its  transit. 

4.1.2.  Ciihtcal  inserts 

The  acoustical  response  of  a  cubical  insert  or  cavity  (sec  Figure  2).  is  similar  to  that  of 
a  rectangular  column.  The  edges  facing  the  incident  acoustic  excitation  are  centres  of  shear, 
(T., ;  see  Figure  7  The  ma  jor  difference  occurs  in  the  o,,.  w  hich  is  of  course  0  for  a  column. 


Table  1 


Materia!  parameters:  hromohuty!  rubber  and  surrounds 


Elastomer 

A 

Steel 

Water 

Young's 

modulus. 

Poisson 

Density. 

1 

Young's 

modulus 

Poisson 

Density 

Density 

/ 

raiic. 

1' 

E 

ratio. 

t> 

t> 

(Pa) 

V 

(kg  m-) 

(Pa) 

\ 

(kg  m') 

(kg/m') 

.S  O  X  10' 

0-4V 

11.70 

1-6  X  10'® 

0-26 

7800 

1000 

Elastomer 

Frequency . 

r 

Relaxation  time  constants 

i') 

Loss  tangent. 

(kHz) 

tan  (>, 

(s  X  10  •') 

(s  X  10  ') 

10 

0-612 

2  i 

6.7 

16 

06.74 

17 

4-.7 

2^ 

06.70 

0-88 

2-6 

40 

0  66V 

()-.76 

1-6 

05 

0  684 

()-.76 

0-99 

100 

069.7 

0-24 

0-6.7 
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l  ipurv  7  Shear  'stress  n,  roiatint:  about  the  i  -axis  around  cubical  cavities,  showing  the  phasing  relationship 
around  and  between  cavities  Lxcitation  Irom  above  figure  in  direction  of  t-axis 


and  which  for  square  cavities  has  maxima  at  the  edges  in  line  with  the  direction  of  advance 
of  the  incident  acoustic  wave.  The  shear  rotating  about  the  .v-axis.  tr...  is  shown  in 
Figure  8. 

4.I..T  Cylindrial  inserts 

The  response  of  a  cylindrical  cavity  excited  normal  to  its  long  axis  is  similar  to  that  of 
a  rectangular  column.  The  shear  stress  pattern  is  shown  in  Figure  9. 


f  igure  S  Shear  stress  o  rotating  about  the  v-axis  around  cubical  cavities,  showing  the  phasing  relationship 
around  and  between  cavities.  Nodal  lines  are  rcadilv  visible.  Excitation  from  above  figure  in  direction  of  .v -axis 
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Figure  9.  Shear  stress  rotating  about  the  r-axis  around  cylindrical  column  cavities,  showing  phasing 
relationship.  Excitation  from  direction  of  I’-axis. 


When  a  cylindrical  air-filled  cavity  is  acoustically  excited  along  its  long  axis,  the  response 
of  the  surrounds  is  similar  to  that  of  a  square  column  for  the  o„  and  a^.  stress  as  expected; 
sec  Figure  10.  The  rotates  about  the  r-axis  and  the  a^.  rotates  about  the  y-axis. 
Additional  motion  occurs  around  the  circumference  of  the  cylinder.  The  associated  shear 
stress  rotating  about  the  .v-axis  is  shown  in  Figure  11.  Conservation  of  momentum 
requirements  require  that  the  basic  response  is  a  two-wavelength  shear  strain  excitation. 
The  presence  of  neighbouring  cavities  is  believed  to  limit  the  cylinder  response  to  this  basic 
mode.  This  result  is  of  importance  in  designing  resonance  absorption  into  acoustic  energy 
absorbing  structures. 


4. 1 .4.  Diamond  and  double  pyramid  cavities 
The  transformation  of  normal  into  shear  stress  is  enhanced  at  the  point  of  these  cavities 
where  the  spatial  gradient  of  the  material  properties  are  greater  than  for  the  cavities 


Figure  ID.  Shear  stress  o,,  rotating  about  the  :-axis  around  the  upper  nm  of  a  cylindrical  cavity,  showing 
the  phasing  relationship.  Excitation  from  above  ligure  m  direction  of  x-axis. 
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Figure  1 1.  Shear  sires<>  a,,  rotating  about  the  .t-axis  around  a  cylindrical  cavity,  showing  phasing  relationship 
The  excitation  was  incident  from  above  figure  along  direction  of  .v-axis. 

considered  so  far.  A  typical  shear  stress  pattern  around  the  point  is  similar  to  that  observed 
at  a  comer  in  square  cavities  and  around  cylindrical  cavities. 

4.1.5.  Solid  inserts 

Typical  shear  stresses  in  and  around  solid,  square,  column  inserts  the  acoustic 
impedance,  but  not  the  speed  of  sound,  of  which  is  equal  to  that  of  the  background  matrix 
material,  are  shown  in  Figure  12.  The  shear  stress  is  zero  at  the  surface  in  accordance  with 
acoustic  theory  and  concentrated  along  the  material  interfaces  where  the  sound  waves  in 


Figure  1 2.  Shear  stress  <r,  rotating  about  the  r-axis  around  an  insert  the  acoustic  impedance,  but  not  the  speed, 
of  which  IS  equal  to  that  of  the  background  matrix  matenal.  Excitation  was  along  the  .v-axis. 
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the  two  media  are  out  of  phase.  Similar  patterns  of  stress  were  computed  for  cubical 
inserts  such  as  shown  in  Figure  2. 


5.  ENERGY  DISSIPATION 

Energy  dissipation  properties  determined  the  efficiency  with  which  an  element  of  a 
given  material  composite  attenuates  sound  and  vibration  and  is  fundamental  for  the 
design  of  sound-attenuating  coatings  and  machinery  mounts.  Irreversible  energy  losses  in 
hysteretic  fluids,  such  as  viscoelastic  elastomers  or  plastic  materials  are  caused  by 
particle-particle  interaction  when  an  element  of  the  fluid  is  distorted  by  the  acoustic  or 
vibrational  stress  wave.  These  losses  are  proportional  to  the  product  of  time  rate  of 
change  of  strain  and  stress,  denoted  by  L  in  equation  (33).  and  hence  are  dominant  at  the 
corners  of  cavities  and  edges  of  inserts  where  both  strain  and  stress  change  rapidly.  This 
general  picture  of  the  energy  dissipation  in  the  material  changes  only  slightly  with 
frequency. 


Table  2 

Relative  power  dissipation  in  composites 

Thickness  of  composite  =  0-5  cm  coating:  insert  shape  is  rectangular  columns:  insert 
size  =  0  5  cm  wide,  spaced  0-5  cm  apart:  facing  material  is  water:  backing  material  is  steel: 
frequency  of  acoustic  excitation  =  100  kHz;  number  of  computed  oscillations  =  3 


Power 

loss, 

Composite.s  L, 

Material  properties 

Background 

Insert 

A 

/ 

tan  6, 

P 

tan^^ 

P 

( - 

/ 

tan  b. 

P 

tan  <5^ 

P 

A 

0-30 

0-854 

0-000 

0-689 

0-010 

0-689 

0-606 

0-0 

0-000 

0-16 

2-0 

B 

000 

0-531 

0-0 

0-000 

0-12 

1-02 

0-527 

00 

0-004 

009 

1-7 

C 

001 

0854 

0-0 

0-689 

0-01 

1-689 

1-334 

00 

0-444 

0-50 

4-0 

D 

2C 

0  854 

00 

0-689 

0-01 

1-689 

0  150 

00 

0-000 

0-16 

2-0 

E 

0  01 

0  444 

00 

0000 

00 

1-0 

1-334 

0-0 

0-444 

0-50 

4-0 

F 

000 

0-444 

00 

0000 

10 

l-O 

1-334 

00 

0-444 

0-50 

1-0 

G 

100 

0824 

0  001 

0209 

0  01 

1  213 

0-604 

0-0 

0-000 

0-16 

2-0 

H 

2300 

0  41*) 

0-002 

0209 

0  01 

1-99 

0  159 

00 

0-000 

0-16 

4-0 

Here  /.  and  pt  are  the  hrsi  and  second  Lame  constants  of  the  material;  tan  i>,  and  tan  are  the  associated 
loss-tangents,  and  p  denotes  the  density  of  the  respective  material.  The  material  properties  are  normalized  with 
respect  to  the  density  and  equivalent  material  modulus  of  water  (I  'p.  where  r  is  the  speed  of  sound  in  water), 
and  the  power  loss  is  normalized  with  respect  to  that  of  material  combination  G.  Thus  all  parameters  in  this 
table  are  dimensionless. 


Composite 

Matenal 

Background 

Insert 

A 

syntactic  foam 

heavv  silicone  rubber 

B 

silicone  rubber 

heavy  silicone  rubber 

C 

syntactic  foam 

heavv  loaded  plastic 

D 

syntactic  foam 

heavv  loaded  silicone  foam 

E 

silicone  rubber 

heavv  loaded  plastic 

I 

silicone  rubber 

plastic 

G 

plastic 

heavv  silicone  rubber 

II 

heavy  plastic 

heavy  loaded  silicone  rubber 
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5.1.  RANKING  OK  ACOUSTIC  POWER  DISSIPATING  COMPOSITES 

Researchers  in  the  acoustic  attenuation  held  have  surveyed  the  range  of  practical,  energy 
dissipative  composites  [1]  which  are  suitable  for  noise  absorbing  structures.  The  composites 
examined  are  shown  in  Table  2.  The  composites  consist  of  elastomer  inserts  in  a  plastic-like 
background  material.  To  illustrate  the  usefulness  of  the  present  model,  it  was  used  rank 
the  average  power  dissipation  properties  of  the  composites  shown  in  Table  2  for  simple 
square  columns.  The  relative  average  power  dissipated  by  the  various  composites  is 
denoted  by  L,  and  is  shown  in  Table  2  [1].  The  average  power  dissipation  L,  is  normalized 
to  the  loss  of  a  composition  believed  to  be  a  practical  composite.  The  calculations 
underlaying  this  ranking  did  not  consider  resonance  effects  due  to  size  or  shape. 
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